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Abstract Distinctive wave forms in the distributions of vertical velocity and temperature of 
both neutral particles and ions are frequently observed from Dynamics Explorer 2 at altitudes 
above 250 km over the polar caps. These are interpreted as being due to internal gravity waves 
propagating in the neutral atmosphere. The disturbances are characterized by vertical velocity 
perturbations of the order of 100 m/s and horizontal wave lengths along the satellite path of about 
500 km. They often extend across the entire polar cap. The associated temperature perturbations 
indicate that the horizontal phase progression is from the nightside to the dayside. Vertical 
displacements are inferred to be of the order of 10 km and the periods to be of the order of 10 3 s. 
The waves must propagate in the neutral atmosphere, but they usually are most clearly 
recognizable in the observations of ion vertical velocity and ion temperature. By combining the 
neutral pressure calculated from the observed neutral concentration and temperature with the 
vertical component of the neutral velocity, an upward energy flux of the order of 0.04 erg/cm 2 - s at 
250 km has been calculated, which is about equal to the maximum total solar ultraviolet heat 
input above that altitude. Upward energy fluxes calculated from observations on orbital passes at 
altitudes from 250 to 560 km indicate relatively little attenuation with altitude. 


Introduction 

The purpose of this paper is to describe some wave motions 
observed over the polar caps at altitudes in excess of 250 km 
using Dynamics Explorer 2 (DE 2) data from the ion drift meter 
(IDM) [Heelis et al., 1981], the retarding potential analyzer 
(RPA) [Hanson ct al.,1981], the wind and temperature 
spectrometer (WATS) [Spencer et a! 1981], and the neutral 
atmosphere composition spectrometer (NACS) [Carignan et 
at-. 1981], The waveforms are interpreted as being due to 
internal gravity waves propagating in the neutral atmosphere. 
The motions are seen as well-formed, quasi -sinusoidal patterns 
of vertical velocity (both ion and neutral) with associated 
temperature and composition perturbations that permit 
evaluation of some properties of the waves; these patterns 
indicate that the horizontal phase progression of the waves is 
from the nightside to the dayside and that the waves transport 
significant amounts of energy upwards, even at altitudes as 
high as the base of the exosphere. Waves of this sort have 
been observed over much of the Earth by Potter et ai. [1976], 
who described their distribution with latitude and possible 
causes. Gross and Huang [1985] have observed waves of 
similar wavelength at low to middle latitudes and altitudes near 
250 km. 
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The observation of internal gravity waves at altitudes as 
high as the exobase was unexpected because of anticipated 
strong attenuation above 120 km due to viscosity and thermal 
conduction. The expectation of strong absorption is based on 
the work of Pittcway and Hines [1963], Midgley and Ucmohn 
[1966], and Hodges [1969]. 

Wave forms are illustrated in Figure 1 for propagating and 
s t a ndin g waves; it will be useful to keep these patterns in mind 
in the subsequent discussion of the observations. The 
displacements from equilibrium positions and the velocities 
associated with the changing displacements are portrayed. In 
the present context the displacements and velocities are the 
components in the vertical direction. From satellites, neither 
the displacements nor the time variations at fixed points in 
space can be observed directly because of the motion of the 
satellite. The satellite observations constitute nearly a 
snapshot at a fixed time because the satellite velocity is so 
high compared to the wave propagation velocity; that is, the 
observations are essentially of spatial variations. The 
vertical displacement perturbations can be inferred from the 
temperature perturbations or from composition changes. The 
displacements and velocity perturbations are in phase 
quadrature in time for both propagating and standing waves. 
For propagating waves the perturbations are also in phase 
quadrature in space, and the displacement profiles lead or lag 
the velocity profiles depending upon whether the coordinate 
axis is in the direction of wave propagation or opposed to it. 
In the case of standing waves the spatial profiles are in phase 
half the time and out of phase half the time, and no 
perturbations occur at the nodes, as illustrated in Figure 1. 
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Figure I. The spatial distribution of displacements and velocities for propagating waves and for standing 
waves at five selected time intervals. At selected instants in time the velocities and displacements for standing 
waves are either in phase (eg., for nfl < \ < x) or in antiphase (eg., for 0 < t < nil) with one another. 


Observations 

Figure 2 shows the wave perturbations seen in the vertical 
ion drifts observed by the IDM on 22 orbital passes over a 
polar cap. Seven-thousand-kilometer (15-min) orbital 
segments are shown, the arbitrarily selected start times being 
indicated in each panel; the time intervals have been chosen 
so that the segments are located over a polar cap. (In three 
cases the initial portion of the segment is missing.) The K p 
magnetic indices are shown in each panel; there is no obvious 
relationship between occurrence of these waves and the K p 
index. However, the interplanetary field was directed 
southward during 33 out of 40 passes on which waves were 
observed and interplanetary field data were available. The 
vertical velocity perturbation amplitudes vary from about 100 
to 300 m/s; other records show smaller amplitudes. Although 
not shown in Figure 2, velocity -correlated ion temperature 
perturbations were present in each case, and these could be 
used to determine whether the horizontal wave progression had 


a component in the direction of satellite motion or opposed to 
it. 

Figure 3 shows IDM data for ion velocity components and 
RPA data for ion temperatures and concentrations from orbit 
8303 on February 1, 1983; the spacecraft altitude was about 
253 km. The oscillations in question extend virtually across 
the entire polar cap. The satellite entered the polar cap near 
local magnetic noon, at 1613:35 UT, and exited near 0300 
MLT, at 1621:09 UT. The waves are seen most clearly in 
terms of the vertical ion velocities, v y . Because the magnetic 
field is almost vertical, ion motions in the horizontal are 
severely constrained and it is no surprise that no evidence of 
wave motions can be seen in the horizontal ion velocity 
components along the orbit or v ? across the orbit- Ion 
temperature measurements are shown in the top panel and the 
oscillations are clearly evident as perturbations of the ion 
temperature. Ion concentrations are shown in the second 
panel, but they show little or no evidence of the oscillations, 
and this is usually the case in other records. 
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Ftoire 2 A selection of records from the ion drift meter on DE 2 showing periodic vertical velocity 
pemubations at altitudes of 250 km and above. The oibit pumber, date* UT at the beginning of the record, 
altitude at the beginning of the record, and K p index are all shown. The ordinate scale is vertical ion velocity 

in meters per second. 


Figure 4 shows the same ion vertical velocity and ion 
temperature records shown In Figure 3. bat for the time 
interval 1614:34 to 1621:14 UT. along with the cross 
correlation between them. The record consists of 400 
observations; the abscissa scale for the middle panel is the 
serial number of the individual observances. The scale for the 
top panel is the displacement to the right of the ion velocity 
H»m relative to the ion temperature data in units of the interval 
between observations. (Positive values are referred to as 
"lead".) To make the cross-correlarioo computation, the 


observations were first convened to perturbation values by 
taking differences from the straight, dotted detrend lines 
shown running through the observations. 

The correlation is greatest for a lag of 11.5 intervals. 
Taking the shifts for the adjoining peaks in the correlation 
function (-65 and +40) as indicative of twice the wavelength, 
the average wavelength along the satellite path is 53.5 
intervals (that is, 53.5 seconds or 414 km). Thus the effective 
phase shift between ion velocity and temperature is 
(11 .5/53.5)360° - 79°. or nearly 90°. Simply overlaying the 
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Figure 3. loo temperatures, concentrations, and velocity components observed from DE 2 on orbit 8303 on 
February l, 1983. The abscissa labels apply to the tic marts at the left-hand ends of the labels. The x axis is 
along the direction of motion of the satellite, the y axis is upwards, and the i axis is perpendicular to the orbit 
plane. Observations are plotted at 1-s (7.74-km) intervals. 


records in the two lower panels also demonstrates dearly that 
the phase difference is about 90*. The exact value of die phase 
shift is not used in the analysis of the vertical component of 
energy flow described below; the phase difference is used only 
to determine whether the horizontal phase progression of the 
waves has a component in the direction of satellite motion or 
opposed to it. 

Because of the predominance of neutrals over ions, the 
waves must propagate in the neutral atmosphere and the 
motions of the ions are driven by the neutral motions. We 
assume that the temperature perturbations have their maximum 
positive values at points of maximum downward displacement 
and thus maximum adiabatic heating. The ion temperature 
perturbations are simply a reflection of the neutral temperature 
perturbations to which they are closely coupled at the altitude 
of observation. Near the points of maximum downward 
displacement, the upward vertical velocities must be near zero 


and increasing with time. The phasing in Figure 4, the upward 
vertical velocity profile lagging the temperature perturbation 
profile, indicates that the horizontal phase progression has a 
component in the direction opposite to that of the satellite 
motion. As the satellite motion in this case was from the 
dayside to the nightside, the direction of horizontal wave 
propagation is from the nightside to the dayside. 

Figure 5 shows the cross correlation between neutral 
vertical velocity and neutral temperature perturbations. In this 
case the observation interval was 8 s, so only 50 observations 
are available during the 400-s period of observation. The 
perturbation values were evaluated in terms of the differences 
from the best-fit dotted detrend lines shown in the figure. The 
cross coneluion is very similar to that for the ion 
temperatures and vertical velocities shown in Figure 4; the 
average phase shift in this case is 95°. 

At tire satellite altitude of about 250 km the ion and neutral 
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Figure 4. Part of tbe data presented in Figure 1 for 400 measurements of the ion vertical velocity and ion 
temperature over the polar cap, and the cross correlation between them* . The dotted lines are (he backgrounds 
from which perturbations were evaluated. The phase shift of nearly 90° is characteristic of propa g a ting waves, 
and the fact that the velocity perturbation lags die temperature perturbation indicates that the horizontal wave 
progression is opposite to tbe direction of satellite motion; that is, propagation is from the nightside to the 
dayside. 


vertical velocities are closely coupled, and there is striking 
observational evidence that this is so. As the m ag ne tic field 
is nearly vertical, the observed ion vertical velocity 
components are essentially the field-aligned components of 
the neutral velocity* The correlation (Le., the cross 

'*( c* .*• 'v 


correlation at zero shift) between the ion and neutral velocity 
perturbations shown in Figures 4 and 5 is 0.93. 

There is good evidence for the presence of horizontal j 
neutral velocity perturbations associated with the waves, tbe 
horizontal perturbation velocities being about one-third as 

* • t. 
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Figure 5. ' Neutral vertical velocities end t emp erat ure s plotted at 8-s intervals over the same time period 
covered by' Figure 4, and the cross correlation between them. The dotted lines indicate the reference levels 
from which the perturbations were evaluated. 
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large as the vertical (e.g., orbits 1182, 1418, 1488, 4698, horizontal winds with strong horizontal gradients, associated 
4828, 6877, 8803). However, in many cases, no related with anti-sunward convection over the polar cap, making it 
perturbations in the horizontal neutral velocities could be necessary to look for small perturbations on steep slopes, 
seen; this might be due to the fact that only the cocnpooent of The correlation (Le., the cross correlation at zero shift) 
the horizontal motion perpendicular to the orbit plane was between the ion and neutral temperature perturbations shown 
observed; the .perturbations oiay^ have been too, small to in Figures 4 and SistOJU. Although the ion and neutral 
recognize. Further; recognition of pert urb a tio n s as so ci a t ed temperature perturbations agree very well, the absolute values 
with the waves, is made difficult by the presence of large are often not in good agreement, probably reflecting 
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inferred from the observed temperature perturbations 
(assuming adiabatic temperature changes due to vertical 
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Upward Energy Transport 
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560 km. In diseases the horiiontal phase propagation of the 
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Figure 8. Five examples of upward energy flow due to waves at altitudes ranging from 310 to 560 km. 
Although the two highest orbits show lesser upward energy flow than those at lower altitudes, the upward 
component of energy flow falls off remarkably little with increasing altitude. 


cases they extended all the way across the polar cap. In a few 
cases there was evidence that the waves extended a short 
distance beyond the polar cap on the dayside. 

A troubling observation is that some examples have been 
found that have clear wave signatures in terms of vertical ion 
velocities but without recognizable accompanying ion 
temperature perturbations. It is a characteristic of standing 
waves that the displacements be zero when the velocity 


perturbations have their maiimum values, so there is some 
temptation to suggest that such observations are indicative of 
standing waves. However this hypothesis seems not to be 
tenable, in that we have found no cases showing clear 
temperature signatures in the absence of clear signatures in 
vertical velocity, something that should be observed with 
equal frequency for standing waves. 

Finally, it is important to note that the closeness of the 
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wave frequency to the Brunt-Vaisala frequency co B , as well as 
the dominance of the vertical component of velocity over' the 
horizontal velocity in the observed wave structures, are natural 
consequences of the dispersion of gravity waves in the upper 
thermosphere. Throughout much of the thermosphere the 
temperature increases with increasing altitude; this causes to B 
to decrease with altitude. In the plane wave approximation, 
Snell's law requires that the horizontal component of the wave 
number be invariant with altitude. Lamb [1945] showed that 
the effective vertical component of the wave number decreases 
with increasing temperature, reaching zero at the level where 
to B - cj , and evanescent beyood. Hence the decay of <o B with 
increasing altitude in the thermosphere imposes a low -pass 
filter with decreasing cutoff frequency on the spectrum of 
gravity wave oscillations. 

The prevalence of observable wave structures near the 
Brunt-Vaisala frequency is apparently owed to the resonance 
condition as <o B approaches , which corresponds to a match of 
the wave frequency with the natural frequency of buoyant 
oscillations of gas parcels. The width of this resonance is 
made finite by viscous dissipation, and gravity waves excited 
by impulsive natural phenomena have a broad spectrum of 
frequencies. As a result, wave-induced atmospheric 
perturbations over a band of wave frequencies near but below 
are selectively amplified by the resonance, creating a 
quasi-sinusoidal structure in which the lower-frequency 
components of the wave spectrum are suppressed. 

In quasi-resonant oscillations near the Bnmt-Vaisala 
frequency the dominance of vertical motion seems obvious. 
An analytic explanation can be found in Figure 3 of Midglcy 
and Licmohn [1966], which shows how the orbits of gas 
parcels in acoustic-gravity waves change with wave 
parameters. What is of interest to the present discussion is 
that the orbits are elliptical, with the major axis nearly 
vertical near the condition co B m to . As to B increases with 
decreasing altitude below about 250 km down to below 100 
km, the major axis of ^the^orbit ellipse /tilts towaid the 
horizontal. Hence the large amplitude oscillations near co B 
that predominate above 250 km are rooted in waves that pass 
through lower altitudes with very classical characteristics, that 
is, with horizontal velocities greater than the vertical. At 
altitudes above about 250, above which the Brunt-Vaisala 
frequency changes relatively little, the vertical velocities are 
naturally greater than the lateral, in agreement with the 
velocity measurements. 
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